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Mott transitions induced by strong electric 
fields are receiving a growing interest. Recent 
theoretical proposals have focused on the Zener 
dielectric breakdown in Mott insulators, however 
experimental studies are still too scarce to con- 
clude about the mechanism. Here we report a 
study of the dielectric breakdown in the narrow 
gap Mott insulators GaTa4Se8-a;Tea^. We find that 
the I-V characteristics and the magnitude of the 
threshold electric field (E^^) do not correspond to 
a Zener breakdown, but rather to an avalanche 
breakdown. 'Eth increases as a power law of the 
Mott Hubbard gap (E^), in surprising agreement 
with the universal law Fith cxE^^^ reported for 
avalanche breakdown in semiconductors. How- 
ever, the delay time for the avalanche that we 
observe in Mott insulators is over three orders 
of magnitude longer than in conventional semi- 
conductors. Our results suggest that the electric 
field induces local insulator-to-metal Mott transi- 
tions that create conductive domains which grow 
to form filamentary paths across the sample. 

In many materials with partially filled d bands, the 
electron motion is frozen as a result of strong repulsive 
interactions, which drive the system into a Mott insulat- 
ing state [U |2]. Introduction of carriers in a Mott in- 
sulator often leads to exotic properties, such as high-Tc 
superconductivity in cuprates [3 and colossal magnetore- 
sistance in manganites [4 . A Mott insulator can also be 
collapsed into a correlated metal by applying pressure. 
Notable examples are V2O3 [5] and organics [6]. An- 
other important type of Mott transition, which is receiv- 
ing a great deal of attention, is the dielectric breakdown 
caused by strong electric fields. Many recent theoretical 
studies have focused on the Zener mechanism of break- 
down for Mott insulators. For instance, calculations were 
performed in ID Hubbard chains using exact diagonal- 
ization [THIQ] and time-dependent density matrix renor- 
malization group [TP , or in the limit of large dimensions 



using dynamical mean field theory [12 . These theoret- 
ical studies predicted a similar non-linear behavior in 
the current- volt age characteristics, and the existence of a 
threshold field (E^^) beyond which a field induced metal 
appears. In the latter case, where the Mott transition 
occurs beyond a critical value of the Coulomb interac- 
tion which is of similar magnitude as the bandwidth, it 
was found that E^^ is proportional to /S? /vp^ with A the 
charge gap, and vp the Fermi velocity. This implies a 
large threshold field of the order of 10^-10^ V/cm, which 
is substantially larger than the experimental values ob- 
served in dielectric breakdown [13 . Surprisingly, com- 
pared to this significant theoretical effort, experimental 
studies of Mott insulators under strong electric fields re- 
main quite scarce. Taguchi et al. observed a dielectric 
breakdown in the quasi-one-dimensional Mott insulators 
Sr2Cu03 and SrCu02 [13 . More recently a compara- 
ble dielectric breakdown was observed in the family of 
chalcogenide Mott insulators AM4Q8 (A = Ga, Ge; M = 
V, Nb, Ta, Mo; Q = S, Se) [T4HT8]. However in both cases 
^th is of the order of 10^-10^ V/cm which does not match 
with theoretical predictions for Zener breakdown. This 
strong discrepancy raises therefore the question of the 
actual physical origin of the observed breakdown. Here 
we focus our study on the chalcogenide Mott insulators 
AM4Q8 and show that the dielectric breakdown is consis- 
tent with an electric field induced avalanche phenomenon 
rather than Zener breakdown. Our work reveals that the 
I(V) curve, the magnitude of ^th and the dependence of 
E^^ versus the Mott-Hubbard gap are those expected for 
an avalanche breakdown, qualitatively similar to what is 
observed in semiconductors. However, our results also 
reveal significant differences, for instance, a delay time 
for the avalanche that is at least three orders of magni- 
tude longer in the Mott insulators AM4Q8 compared to 
the case of conventional band insulators. 

RESULTS 

Mott-Hubbard gap evolution in the 



GaTa4Se8-a;Tea^ compounds 



The lacunar spinel compounds AM4Q8 , with tetrahedral 
transition metal clusters (see Fig. [l]) [19 , represents an 
interesting family of narrow gap Mott insulators. The low 
Mott gap value makes these compounds highly sensitive 
to doping or pressure. Doping of GaV4S8 by Ti leads, 
for example, to a filling control Insulator to Metal Tran- 
sition (IMT) with the emergence of a ferromagnetic half 
metal ^204i21j. On the other hand, an external pressure 
induces a bandwidth controlled IMT and superconduc- 
tivity in GaTa4Se8 [22l|23]. Recently, we have reported 
that electric field leads to a striking volatile IMT or di- 
electric breakdown in GaTa4Se8 [14-17 which is qualita- 
tively different from the non-volatile resistive switching 
based on electrochemical effects reported so far in many 
transition metal oxides p!7| l24| [25] . 

In order to study the origin of this dielectric break- 
down and to explore its relationship with the Mott Hub- 
bard gap we have investigated the electric field effects in 
the series of compounds GaTa4Se8-a?Tea; (x=0, 4.0, 6.0, 
6.5). These compounds were synthesized as described 
elsewhere [26 . Optical conductivity and resistivity mea- 
surements were subsequently undertaken to evaluate the 
gap evolution in this series. The room temperature opti- 
cal conductivity of GaTa4Se8, shown in Fig. [TJd, displays 
all the expected features of a narrow gap Mott Insula- 
tor with a Hubbard energy (U) of the order of 0.55 eV 
[27] . The low bias resistivities p(T) measured with a stan- 
dard four probes technique for the series of compounds 
GaTa4Se8-xTea; exhibit a semiconducting-like behavior 
with a temperature dependence consistent with an acti- 
vation law at high temperature (see Fig. Il}c). The sys- 
tematic change of the activation energy observed in Fig. 
[1] suggests that the substitution of Se per Te permits the 
bandwidth-control of the Mott Hubbard gap (from 100 
to 300 meV) within the GaTa4Se8-a:Tea; series of com- 
pounds, thus enabling the study of the threshold field 
dependence with the size of the gap. 

Dielectric breakdown experiments in the 

GaTa4Se8-a:Tea; Mott insulators 

Dielectric breakdown was systematically studied in 
this series of narrow gap Mott insulators using a pro- 
cedure described in the Methods section. Figure 2 shows 
that these compounds undergo a sudden volatile dielec- 
tric breakdown above a certain threshold electric field 
(^1-10 kV/cm). As an example. Fig. 2c depicts the 
typical temporal evolution of the voltage V5'(t) across a 
GaTa4Se8 crystal while a series of short (10 to 70 jas) volt- 
age pulses {V pulse) are applied. When the electric field 
applied to the crystal exceeds a threshold field of ^ 2.5 
kV/cm the crystal systematically undergoes a fast tran- 
sition to a low resistance state (i.e. from 70 kQ to 5 kQ ) 
called a dielectric breakdown or resistive switching. This 
transition is volatile and perfectly reproducible. In fact. 
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FIG. 1: Resistivity and optical conductivity measure- 
ments, (a) representation of the lacunar spinel structure 
of the AM4Q8 compounds. Blue, orange and Yellow filled 
circle represents M, A and Q atoms respectively. The M4 
tetrahedral cluster and the AQ4 tetrahedra are highlighted in 
blue and orange, respectively (b) experimental optical con- 
ductivity of GaTa4Se8 (green line) showing that the optical 
conductivity is that expected for a Mott Insulator with an 
Hubbard energy U= 0.55eV and a gap close to 0.12 eV. The 
gray solid line corresponds to the fit of optical conductivity 
spectrum using two contributions : a broad mid-infrared band 
(red shaded area) and high-frequencies contributions (thick 
solid red line), (c) Resistivity as a function of inverse tem- 
perature for GaTa4Se8 (filled green triangle), GaTa4Se4Te4 
(black open triangle), GaTa4Se2Te6 (blue open square), and 
GaTa4Sei.5Te6.5 (filled red square). 



after the pulse application the low-bias resistance returns 
to its original value. The delay t delay between the appli- 
cation of the pulse and the dielectric breakdown shifts 
to smaller values with increasing applied voltage. The 
voltage drop on the sample V^ in the breakdown state 
is always Vs ^ 9 V (or Eth ~ 2.5 kV/cm) whatever the 
voltage applied to the sample before the resistive switch- 
ing. On the other hand, the current flowing through the 
sample after the dielectric breakdown shows a highly non- 
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FIG. 2: Resistive switching experiments and model- 
ing. (a)I-V curves extracted from pulse (open circles and 
squares) experiments for a crystal of GaTa4Sei.5Te6.5. The 
circuit used to apply voltage pulses is shown in inset, (b) 
I-V curves extracted from pulse (open circles and squares) 
or DC experiments (plain or dashed lines) for a crystal of 
GaTa4Se8. (c) Evolution of the voltage Vs(t) across a crystal 
of GaTa4Se8 when short voltage pulses Ypulse (10-70 /xs) of 
increasing magnitude are applied to the circuit. For example, 
Ypulse = 5V , 35V, 40V, 44V, 50V, GOV and 86V for the 
gray, blue, black, green, pink, light blue and orange curves, 
respectively. These curves show a clear evolution of the time 
{t delay) at which occurs the transition to a low resistance state 
(resistive switching) as a function of Y pulse- In (a) and (b) 
the I-V curves are extracted from Vs(t) curves similar to the 
one shown in (c) and from the corresponding I(t) curves not 
shown here. In (a) and (b) open circles correspond to the 
I-V values measured before the resistive switching and open 
squares correspond to the I-V values measured after the re- 
sistive switching, (d) Results of mathematical modeling of 
Vs(t) using a one dimensional resistive network. The time at 
which occurs the resistive switching decreases as the applied 
voltage to the resistor network is increased from the black to 
the orange curve. 



ohmic behavior in the transited state in both, pulse or dc 
measurements (see Fig. 2-b). Thus, we infer that Eth (i.e. 
Yth/d where d is the distance between electrodes) corre- 
sponds to the breakdown field for the limiting case tdeiay 
-^ oo. The same qualitative behavior was observed for 
GaTa4Sei.5Te6.5 (see Fig. 2a), all other GaTa4Se8-xTea^ 
compounds, and more generally for all AM4Q8 (A = Ga, 
Ge; M = V, Nb, Ta, Mo; Q = S, Se) Mott insulator 
compounds [17j. 



Dielectric breakdown mechanism in the 

GaTa4Se8-a:Te2^ Mott insulators 

For conventional semiconductors the dielectric break- 
down is well documented and different mechanisms are 
proposed such as Zener tunneling [28] or avalanche pro- 
cess [29 . The threshold field observed in the AM4Q8 
compounds is at least 3 orders of magnitude smaller than 
the threshold field of the order of 10^ - lO^V/cm expected 
for Zener breakdown in Mott insulators [7 . This defini- 
tively rules out a Zener breakdown mechanism. Con- 
versely, the magnitude of the threshold field in AM4Q8 
Mott insulators compares well with the threshold field 
values of 1 and 40 kV/cm observed for avalanche break- 
downs in narrow gap semiconductors InSb and In As, re- 
spectively [30l|3T]. The I(V) characteristics observed in 
the AM4Q8 compounds and in semiconductors like InSb, 
InAs, Si or Ge are also very similar [29]. Indeed, all 
these compounds exhibit a steep increase of the current 
at the threshold field, which is almost vertical. In narrow 
gap semiconductors this is related to the impact- ion- 
ization process. An avalanche breakdown indeed occurs 
when the electrons accelerated by the applied voltage 
can gain enough energy to ionize by impact other elec- 
trons above the band gap and generate electron-hole pairs 
[29j . When the creation of electron- hole pair by impact- 
ionization overcompensates the electron-hole recombina- 
tion process (Auger Thermal process) a multiplication of 
carriers occurs and induces a dielectric breakdown with 
a steep increase of the current at the threshold field. For 
high current density the avalanche breakdown is charac- 
terised by a negative differential resistance (NDR) in the 
current-voltage characteristic which is the consequence 
of the contribution of holes to the field distribution along 
the sample \^ . As shown in Fig. 2, such a NDR effect 
is observed on the I-V curve of GaTa4Se8 measured in 
DC mode using a current source. Again, this strongly 
supports an avalanche mechanism as the origin of the 
dielectric breakdown. 

Dependence of the avalanche threshold field in 
Mott insulators and classical semiconductors 

For semiconductors the avalanche threshold field varies 
as a function of a power law of the band gap [29] and 
follows the universal law Efh oc Eg^-^ [30[[3T]. For com- 
parison Fig. 3 displays the variation of the threshold 
field as a function of the Mott-Hubbard gap for all the 
AM4Q8 compounds. This figure reveals a remarkable 
power law dependence of the threshold fields with the gap 
E^^ (X Eg^-^ (see blue curve in fig 3). The inset of Fig. 3 
compares also the threshold fields versus gap dependence 
for the AM4Q8 compounds and classical semiconductors. 
This comparison reveals that for all Mott and band in- 
sulators or semiconductors the threshold fields line up 
along a universal law E^^ ex Eg^*^. This is again strong 
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FIG. 3: Dependence of Eith in Mott insulators 

and semiconductors. Threshold Electric field (inducing 
avalanche breakdown) as a function of the Mott gap for vari- 
ous AM4Q8 compounds. The solid blue curve correspond to 
a power law dependence Eth oc Eg^"^. Inset : comparison 
of the threshold fields versus gap dependence for the AM4Q8 
compounds and for classical semiconductors. The solid blue 
line displays the universal law Et/i[kV/cm] = 173 (E^[eV])^"^ 
observed for semiconductors. 



evidence that the breakdown observed in the Mott Insu- 
lators AM4Q8 originates from an avalanche phenomenon. 

DISCUSSION 

This scaling of the Mott- and band-insulators along 
the same universal law E^/i = 1.73 E^^-^ indicates that 
the threshold field in the avalanche breakdown depends 
only on the gap whatever its origin and nature. This 
is rather surprising as the excitation gap in Mott insu- 
lators arises from electron-electron repulsions, (i.e. from 
strong correlation effects) , in stark contrast with the case 
of conventional semiconductors. Despite this similarity, 
as we mentioned before, the main qualitative difference 
between the electric breakdown in semiconductors and 
Mott insulators is the significantly longer tdeiay of the 
latter. The origin of the difference may be due to the 
different coupling to the lattice degrees of freedom. The 
Mott transition is indeed accompanied of a compress- 
ibility anomaly [32l [33] which leads to a large volume 
change [34]. In fact, the giant electromechanical cou- 
pling observed during an STM studies by applying elec- 
tric pulses between the tip and a cleaved GaTa4Se8 crys- 
tal surface supports this scenario [16]. Therefore, the 
strong coupling of the electronic and lattice degrees of 



10° - 



GaNb4Se8^^^ ^ 



GaTa4Sei5Tee5 

GaTa4Seg 

GaV,S« 



10° 



10^ 10^ 

E-Eth (kV / cm) - V 




V, (a.u.) 



FIG. 4: Avalanche velocity in Mott Insulators, (a) 

Avalanche velocity v = d/ tdeiay for various AM4Q8 com- 
pounds as a function of excess electric field E-Et/^ where E 
is the electric field applied on the sample, (b) velocity of the 
transition 1 /tdeiay extracted from the mathematical modeling 
of the resistive switching shown in figure 2d. 



freedom is likely to slow down significantly the velocity 
of the avalanche propagation. Interestingly, Fig. 2 shows 
that the dielectric breakdown in the AM4Q8 Mott insu- 
lators occurs with a ideiay in the jis range (from a few 
/is to the ms) for typical sample lengths (i=20-40 jiia. In 
contrast, in semiconductors the tdeiay lies in the ps to the 
ns range, that is, several orders of magnitude faster J29] . 
Figure 4 displays the avalanche velocity v = d/tdeiay as 
a function of excess electric field E-E^/^ where E is the 
electric field applied on the sample. This figure shows 
that V varies as (E-Et^)^. This remarkable power law de- 
pendence as well as the magnitude of the velocity in m/ s 
range recall a phase growth mechanism [35| or the prop- 
agation of dislocations [36^ 37 . This analysis therefore 
suggests that the excess electric field triggers the local 
growth of conducting domains where a Mott transition 
has taken place, which eventually form filamentary paths 
across the sample as the avalanche breakdown occurs. 

A mathematical modeling of the previous scenario was 
undertaken to describe the experimental data. This 
model considers a one dimensional resistive network 
where each cell may be in one of two states, insulator 
(I) or metal (M), characterised by their corresponding 



resistance r/ and tm, and their energy Ej and Em- This 
model also assumes that the two states are separated by 
an energy barrier of height Eb-Ej. The energy landscape 
used for the calculation is shown in Supplementary Fig- 
ure SI. Initially, all cells are in the insulating state, since 
the system is a Mott insulator. An external voltage pulse 
Vpuise is then applied on the resistor network and the lo- 
cal voltage drops on the cells are computed. The cells 
may undergo a transition from state I to M by a thermal 
activation process. The effective activation energy barrier, 
for an I to M transition of the i-th cell under the ac- 
tion of the applied voltage, is(EB - E/ - AV^), where 
AVi is the computed voltage drop at the cell. These as- 
sumptions lead to non-linear differential equation for the 
rate of cell transitions per unit time. Further details of 
the mathematical expressions are provided in the Meth- 
ods section. The solution of the differential equation is 
easily obtained by standard numerical integration. The 
results are shown in Fig. 2(d). We observe that this sim- 
ple model qualitatively captures the behavior of tdeiay of 
the avalanche mechanism. Initially, the conductive path 
modeled by the resistive network is highly resistive, so 
there is a small voltage drop on each cell, and the rate of 
transitions is low. Successive transitions, however, lower 
the resistance of the path, as tm « r/, therefore in- 
creasing the current and also the voltage drop AV^ in 
the cells that still remain insulating. The increased AV^ 
lowers the local energy activation barriers, thus augment- 
ing the transition rates. A larger fraction of cells undergo 
a transition and, eventually, a runaway process take place 
were the rest of the cells rapidly completes their transi- 
tions. Interestingly, this simple model also qualitatively 
accounts for the field dependence of the velocity of the 
transition as shown in Supplementary Figure S2 and in 
Fig. 4. However the power law behavior predicted by 
this simple model is parameter dependent, i.e. is non 
universal. 

To conclude we have reported for the first time an 
avalanche breakdown in the narrow gap Mott insula- 
tors AM4Q8. The phenomenology of the breakdown 
(I(V) characteristic) and the magnitude of the thresh- 
old electric field are comparable with those reported for 
avalanche breakdowns in semiconductors. Moreover the 
threshold field dependence with the width of the Mott 
Hubbard gap shows a surprising agreement with the uni- 
versal law Efh = 1.73 Eg^-^ reported for semiconductors, 
despite the fact that the AM4Q8 are strongly correlated 
systems. Our experimental and theoretical investigations 
suggest that, unlike semiconductors, avalanche break- 
down in Mott insulators is associated to the growth of 
conducting filamentary paths of a metallic metastable 
phase within the sample. The strong coupling of elec- 
tronic and lattice degrees of freedom due to compress- 
ibility anomalies may account for the reduced avalanche 
speed observed in the Mott insulators compared to con- 
ventional band insulators. Finally, our work suggests 



that avalanche breakdown may be a general feature of 
narrow gap Mott insulators. 

METHODS 

Sample preparation 

Single crystals of GaTdi^Ses-yTey with y ranging from 
to 6.5 were prepared using stoichiometric mixtures of 
elemental gallium, tantalum, selenium and tellurium (pu- 
rities > 99,5 %) [26 . These mixtures were loaded in evac- 
uated sealed quartz ampules of about 10 cm and heated 
up at 300° C/h to 1000° C and hold at this temperature 
during 160h. The quartz tubes were subsequently cooled 
down, first slowly to 800° C at l°C/h and then faster 
(300° C/h) to room temperature. These syntheses yielded 
black powders containing a high yield of GaTdi^Ses-yTey 
metallic-gray cubic and tetrahedral crystals (typical size 
< 300 /im) along with a small amount of impurity phases 
(e.g. GaTe, TaSe2). Large amounts of tetrahedral or 
cubic-shaped crystals of all other members of the AM4Q8 
family have been obtained using a similar procedure. 

Transport measurements 

The AM4Q8 crystals used for transport measurements 
were contacted using 17 /im gold wires and carbon paste 
(Electrodag PR-406), and then annealed in vacuum at 
150°C during 30 minutes. The low-bias resistance of the 
AM4Q8 was measured using a source- measure unit (ei- 
ther a Keithley 236 or a Keithley 6430) by a standard 
two- or four-probe technique. We checked that the con- 
tact resistances were much smaller than the sample re- 
sistance. Voltage pulses were applied using an Agilent 
8114A. During the pulse, the voltage and current across 
the sample were measured with a Tektronix DPO3034 
oscilloscope associated with a IeS-ISSD210 differential 
probe, using the simple circuit described in the inset of 
Fig 2-a. 

Mathematical Modeling 

In order to obtain Figures 2d and 4b we numerically 
integrated 



dx 

—-= ^ - exp 
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(1) 



by using 4th order Runge-Kutta method (time step = 
0.001). Here x is the fraction of cells in metallic state. 
Eb, E/, and AV^ are defined as displayed in Supplemen- 
tary Figure SI. We set v = 0.001, ksT = 1,Eb-Ei = 1 
and 



AVi = K. 



ri 



ri - (I - x) ^ vm - X ^ Rl 



(2) 



For Figure 2d we simulated Vpuise = 2.5, 3, 3.5, 4, 4.5, 5 
and 5.5. Figure 4b presents the switching times for the 5 



higher values. The ID resistivities were set as 10 • tm = 
Rl and r/ = 10 • Rl- Vs was obtained as: 



Vs = Vp, 



ri - {1-x) ^tm -x 
ri - (I- x) ^vm - x^ Rl 



(3) 



t 



delay 



was calculated as the inverse of the number of 



Runge-Kutta cycles needed for x to arrive to 1. An- 
other output of this numerical simulation is presented in 
Supplementary Figure S2, the evolution of ^ in time. 

ACKNOWLEDGMENTS 

This work was supported by the French Agence Na- 
tionale de la Recherche through the funding of the NV- 
CER (ANR-05-JCJC-0123-01) and NanoMott (ANR-09- 
Blan-0154-01) projects. 

Author Contributions 

L.C., E.J. and B.C. initiated the project. V.G., L.C., 
E.J. and B.C. performed the synthesis and the carac- 
terisations of the crystals, and the resistive switching 
experiments. V.T.P. did the optical conductivity mea- 
surements. P.S. and M.R. did the modeling. L.C., E.J., 
B.C., V.G., V.T.P., D.R., T.C., P.S. and M.R discussed 
and wrote the article in the framework of the NanoMott 
project. 

Competing financial interest 

the author declare no competing financial interest 
REFERENCES 



Laurent.Cario@cnrs-imn.fr 

Mott, N.F. Metal-Insulator Transitions. Taylor & Fran- 
cis, London (1990) 

Imada, M. et al. Metal-insulator transitions, Rev. Mod. 
Phys. 70, 1039-1263 (1998). 

Bednorz, J. G. et al. Possible high Tc superconductiv- 
ity in the Ba-La-Cu-0 system, Z. Physik B 64, 189-193 
(1986). 

Jin, S. et al. Thousandfold Change in Resistivity in Mag- 
netoresistive La-Ca-Mn-0 Films, Science 264, 413-415 
(1994). 

Limelette, P. et al. Universality and Critical Behavior at 
the Mott Transition, Science 302, 89-92 (2003). 
Kagawa, F. et al. Unconventional critical behaviour in 
a quasi- two-dimensional organic conductor. Nature 436, 
534-537 (2005). 

Oka, T. et al. Ground-State Decay Rate for the Zener 
Breakdown in Band and Mott Insulators, Phys. Rev. Lett. 
95, 137601 (2005). 

Oka, T. et al. Dielectric breakdown in a Mott insulator: 
Many-body Schwinger-Landau-Zener mechanism studied 
with a generalized Bethe ansatz, Phys. Rev. B 81, 033103 
(2010). 



[9 

[lo: 
[11 

[12: 
[13: 

[14 

[15: 

[16 

[17 

[is: 

[19 

[20: 
[21 
[22: 

[23: 

[24 
[25: 

[26: 

[27 

[28: 
[29: 



Oka, T. et al. Breakdown of a Mott Insulator: A Nonadi- 
abatic Tunneling Mechanism, Phys. Rev. Lett. 91, 066406 
(2003). 

Oka, T. et al. Nonlinear transport in a one-dimensional 
Mott insulator in strong electric fields. Physica B: Cond. 
Matter 359-361, 759-761 (2005). 

Heidrich-Meisner, F. et al. Nonequilibrium electronic 
transport in a one-dimensional Mott insulator, Phys. 
Rev. B 82, 205110 (2010). 

Eckstein, M. et al. Dielectric Breakdown of Mott Insula- 
tors in Dynamical Mean-Field Theory, Phys. Rev. Lett. 
105, 146404 (2010). 

Taguchi, Y. et al. Dielectric breakdown of one- 
dimensional Mott insulators Sr2Cu03 and SrCu02, Phys. 
Rev. B 62, 7015 (2000). 

Vaju, C. et al. Electric-pulse-driven electronic phase sep- 
aration, insulator-metal transition, and possible super- 
conductivity in a Mott insulator. Adv. Mater. 20, 2760- 
2765 (2008). 

Vaju, C. et al. Electric-Pulse-Induced Resistive Switch- 
ing and Possible Superconductivity in the Mott Insulator 
GaTa4Se8. Microelectronics Eng. 85, 2430-2433 (2008). 
Dubost, V. et al. Electric-Field- Assisted Nanostructuring 
of a Mott Insulator. Adv. Funct. Mater. 19, 2800-2804 
(2009). 

Cario, L. et al. Electric-Field-Induced Resistive Switch- 
ing in a Family of Mott Insulators: Towards a New Class 
of RRAM Memories. Adv. Mater. 22, 5193-5197 (2010). 
Souchier, E. et al. First evidence of resistive switching in 
polycrystalline GaV4S8 thin layers. Phys. Status Solidi 
RRL6, 53-55 (2011). 

Ben Yaich, H. et al. Nouveaux chalcogenures et chalco- 
halogenures a clusters tetraedriques Nb4 ou Ta4 , J. Less- 
Common Met. 102, 9-22 (1984). 

Vaju, C. et al. Metal-metal bonding and correlated 
metallic behavior in the new deficient spinel Gao.87Ti4S8. 
Chem. Mater. 20, 2382-2387 (2008). 
Dorolti, E. et al. Half-Metallic Ferromagnetism and Large 
Negative Magnetoresistance in the New Lacunar Spinel 
GaTi3VS8. J- Amer. Chem. Soc. 132, 5704-5710 (2010). 
Abd-Elmeguid, M.M. et al., Transition from Mott Insula- 
tor to Superconductor in GaNb4Se8 and GaTa4Se8 under 
High Pressure, Phys. Rev. Lett. 93, 126403 (2004). 
Pocha, R. et al. Crystal Structures, Electronic Prop- 
erties, and Pressure-Induced Superconductivity of the 
Tetrahedral Cluster Compounds GaNb4S8, GaNb4Se8, 
and GaTa4Se8, J. Am. Chem. Soc. 127, 8732-8740 
(2005). 

Waser, R. et al. Nanoionics-based resistive switching 
memories, Nature Mater. 6, 833-840 (2007); 
Rozenberg M.J. Resistive switching. Scholarpedia 6(4) 
11414 (2011). 

Guiot, V. et al. Control of the electronic properties and 
resistive switching in the new series of Mott Insulators 
GaTa4Se8-yTe^ (0<y<6.5). Chem. Mater. 23, 2611-2618 
(2011). 

Ta Phuoc, V. et al. Optical Conductivity Measure- 
ments of GaTa4Se8 Under High Pressure: Evidence of 
a Bandwidth- Controlled Insulator-to-Metal Mott Tran- 
sition. Phys. Rev. Lett. 110, 037401 (2013). 
Zener, C. , A Theory of the Electrical Breakdown of Solid 
Dielectrics, Proc. Royal Soc. A 145, 523-529 (1934). 
Levinshtein, M. et al. Breakdown Phenomena in Semi- 
conductors and Semiconductor Devices, World Scientific 



Publishing Co. Pte. Ltd., Singapore, (2005). 

[30] Hudgins, J. L. et al. An Assessment of Wide Bandgap 
Semiconductors for Power Devices, IEEE Transactions 
on Power Electronics 18, 907-914 (2003). 

[31] Hudgins, J. Wide and Narrow Bandgap Semiconductors 
for Power Electronics: A New Valuation, J. Electronic 
Mater. 32, 471-477 (2003). 

[32] Kotliar, G. et al Compressibility Divergence and the Fi- 
nite Temperature Mott Transition, Phys. Rev. Lett. 89, 
046401 (2002). 

[33] Hassan, S. R. et al. Sound Velocity Anomaly at the 
Mott Transition: Application to Organic Conductors and 



V2O3, Phys. Rev. Lett. 94, 036402 (2005). 
[34] McWhan, D.B. et al. Metal-Insulator Transitions in Pure 

and Doped V2O3, Phys. Rev. B 7, 1920-1931 (1973) 
[35] Willnecker, R. et al. Dislocation velocities, dislocation 

densities, and plastic flow in lithium fluoride crystals. 

Appl. Phys. Lett. 56, 324-340 (1990). 
[36] Messerschmidt, U. Dislocation Dynamics During Plastic 

Deformation, Springer- verlag, Berlin (2010). 
[37] Johnston, W. G. et al. Dislocation velocities, dislocation 

densities, and plastic flow in lithium fluoride crystals, J. 

Appl. Phys. 30, 129-145 (1959). 



